
mata:
	pointer vector callf_setup(pointer(matrix function) f,real scalar n,a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13,a14,a15)
	{
		pointer vector p
		real scalar mycopy
		p=J(1,15+2,f)
		mycopy=n;p[2]=&mycopy
		p[3]=&a1;p[4]=&a2;p[5]=&a3;p[6]=&a4;p[7]=&a5;p[8]=&a6;p[9]=&a7;p[10]=&a8;p[11]=&a9;p[12]=&a10;p[13]=&a11;p[14]=&a12;p[15]=&a13;p[16]=&a14;p[17]=&a15
	
		return(p)
	}
end
